Pinning of dynamic spin density wave fluctuations 
in the cuprate superconductors 
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We present a theory of the pinning of dynamic spin density wave (SDW) fluctuations in a d-wave 
superconductor by local imperfections which preserve spin-rotation invariance, such as impurities or 
vortex cores. The pinning leads to static spatial modulations in spin-singlet observables, while the 
SDW correlations remain dynamic: these are the 'Friedel oscillations' of a spin-gap antiferromagnet. 
We connect the spectrum of these modulations as observed by scanning tunnelling microscopy to 
the dynamic spin structure factor measured by inelastic neutron scattering. 



I 



I 

o 
o 



> 

o 

O 



I 

O 

o 



X 



Many studies of the cuprates assume that high temper- 
ature superconductivity is best understood by a theory of 
doping mobile charge carriers into a Mott insulator which 
is paramagnetic «-e. a Mott insulator with no static 
spin moment on any site and dynamic antiferromagnetic 
spin density wave (SDW) correlations. Studies of param- 
agnetic Mott insulators on the square lattice of Cu ions 
1^ showed that the most promising candidates have bro- 
ken translational symmetry (and confinement of = 1/2 
spinon excitations) associated with the appearance of 
spontaneous bond charge order. Here, and henceforth, 
we define the term "charge order" very generally: there 
is a periodic spatial modulation in observables which are 
invariant under spin rotations and time reversal, such as 
local electron/hole density of states (LDOS), the spin ex- 
change energy or the electron kinetic energy contained in 
a link of the square lattice; the modulation in the total 
charge density may well be unobservably small because 
of screening by the long-range Coulomb interactions. 

Upon doping the charge-ordered paramagnetic Mott 
insulator |^ , superconductivity co-exists with charge or- 
der for a finite range of carrier concentrations (S), and a 
d-wave superconductor with the full square lattice sym- 
metry appears above a critical S. Almost all theoretical 
and experimental studies have concluded that the low 
temperature properties of this d-wave superconductor at 
optimal S are qualitatively identical to those of a BCS 
superconductor obtained by the Cooper instability of a 
metallic Fermi liquid. This conclusion raises the ques- 
tion of whether the connection to the Mott insulator is 
even necessary in a description of the low temperature 
properties of the superconductor at or above optimal S. 

A number of works |^ argued that memory of the Mott 
insulator should survive in and around vortices in super- 
conducting order: the suppression of superconductivity 
in the vortex core implies that the Cooper pairs are not 
condensed, but the electrons should still strive to retain 
the exchange correlation energy of the Mott insulator. 
This reasoning, and the studies of the doped paramag- 
netic Mott insulator 0, led to the suggestion 0] that 



the static charge order, along with dynamic SDW corre- 
lations, should appear near vortex cores. 

Hoffman et al. ^ introduced a novel scanning tun- 
nelling microscopy (STM) technology of atomically reg- 
istered spectroscopic mapping. Applied to slightly over- 
doped Bi2Sr2CaCu208+2;, they detected LDOS modu- 
lations around the vortex cores at wavevectors — 
(7r/2, 0) and 'Key — (0,7r/2) (the square lattice spac- 
ing has unit length and this modulation has a period 
of four lattice spacings), co-existing with well estab- 
lished superconductivity Bulk charge order with 
this, and related, periods has been discussed in insulat- 
ing/superconducting paired hole states from a number 
of different theoretical perspectives [|[ |, H ||, |l^, |l|, 
p^ . Neutron scattering studies of the optimally doped 
cuprates have not (so far) seen dynamic or static charge 
order, but do see coUinear SDW fluctuations at wavevec- 
tors — (37r/4,7r) and K^j^ — (tt, 37r/4) and at en- 
ergies above a spin gap A. On symmetry grounds [^ , 
collinear SDW correlations are accompanied by charge 
order correlations at wavevectors 2'K.ax,y, which equal the 
values of itKcaj.j,, modulo reciprocal lattice vectors. This 
connection was used to describe the spatial envelope 
of the charge order in the doped Mott paramagnet by the 
pinning of a 'sliding' degree of freedom of the SDW (the 
spin-rotation degree of freedom is not pinned and remains 
dynamic) by imperfections, such as vortex cores. 

This paper presents a description of the energy depen- 
dence of the LDOS measured in STM in a simple Gaus- 
sian model of the pinning of dynamic SDW fluctuations. 
We connect the spectroscopic information contained in 
the STM observation at wavevectors Jicx,y to the inelas- 
tic neutron scattering spectra at wavevectors K.sx,y Our 
theory applies to pinning by impurities and most other 
imperfections that preserve spin rotation invariance. It 
also applies to pinning by vortex cores, but neglects the 
Doppler shift in the quasiparticle energies induced by 
the superflow ||l^: this shift has no contribution at the 
wavevectors 'Kcx,y, and so its neglect is justifiable. Our 
results should be distinguished from other recent propos- 



als [|16| which assume static SDW order, for which there 
is no direct evidence in the system of Ref . ||. 

We begin by identifying the complex order parameters, 
= x,y,z extends over the directions in spin 
space) for SDW fluctuations; the spin operator on the 
lattice site r is 

S^{r,T) = Re [e'K- (r, r) +e^K.,.r^^^(j.^^)] 

where r is imaginary time. We use a Gaussian action for 
^x,ya fluctuations: 

Si^T J2 X^'(q,^n)|$.a(q,w„)|' + (x^y); (1) 

we have Fourier transformed from (r, r) to wavevectors 
q and Matsubara frequencies LOn at a temperature T. 
The dynamic spin susceptibilities Xx,y will be inputs from 
which we will deduce the dynamic structure of the pinned 
charge order. However, this should not be interpreted as 
a causal connection from the SDW fluctuations to the 
charge order; rather, our method relies on the assumed 
proximity of the superconductor to a quantum transition 
to a superconducting state with long-range SDW order 

1^. The presence of such a quantum critical point, and 
the values of 'K.sx,y near it, are determined by a compli- 
cated interplay between spin and charge correlations in 
the superconductor and the proximate Mott insulator. A 
corresponding theory using the charge density wave order 
parameter can also be developed near a charge-ordering 
transition in the superconductor; this is not discussed 
here and it yields results closely related to those in [ p^ . 
We will take a simple functional form for Xx,y (see be- 
low) , with a spectral weight which vanishes below a spin 
gap energy A (a more accurate form can, in principle, 
be extracted from neutron scattering experiments). The 
value of A can be tuned by an applied magnetic field per- 
pendicular to the square lattice layers , and its primary 
effect is a reduction in A. 

Invariance under translation by Ux lattice spacings in 
the X direction (for arbitrary integer n^) leads to the 
symmetry ^x,a — > e^'^^^^'^^xa obeyed by 5$ (and simi- 
larly for ^ya)- As long as this symmetry is preserved, no 
observable will exhibit static charge order at wavevectors 
Kcx.y However, this sliding translation symmetry is bro- 
ken by any imperfection (vortex core or impurity) that 
may be centered near r — Tq, provided the spatial extent 
of the imperfection is not much larger than the charge 
order period (4 lattice spacings). We account for such an 
imperfection by adding a pinning term to the action 

5pin = -Y.h^ [C.*L(ro, T) + {x~*y) + c.c] , (2) 

where (^x,y are complex coupling constants representing 
the pinning potential. Note that 5pin contains the sim- 
plest terms which break the ^x.a e^"-^'"^^^xa sym- 
metry, but preserve spin rotation invariance. Additional 
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FIG. 1; One-loop diagrams for the self energy of the host 
conduction electrons. Thick lines are bare conduction elec- 
tron propagators, a single/double thin line denotes a bare/full 
SDW propagator, and the cross represents the scattering of 
the SDW fluctuations by the imperfection pinning term. 

terms like ^xaiT^o, T)$j,a(ro, r) are also permitted but we 
will not include them for simplicity; such terms lead to a 
90° turn in an SDW fluctuation and can be expected to 
have a small amplitude on physical grounds — their pres- 
ence will lead to additional charge order at iKsx ± 
We also note that the theory for the superconductor prox- 
imate to a charge ordering transition has a pinning term 
linear in the charge order parameter. 

STM involves tunneling of single electrons, and so to 
compute the STM spectra we couple the Gaussian SDW 
fluctuations described by iS$-|-iSpi„ to the electrons, which 
we describe by a standard d-wave BCS model: 

Hbcs = E *k [(^k - /i)r^ + AkT^] *k. (3) 

k 

Here ~ (ckticLkj) is a Nambu spinor at momentum 
k = (kx,ky), r^'J*'^ are Pauli matrices in particle-hole 
space, and is a chemical potential. For the kinetic 
energy, Ck we have first (t) and second (f) neighbor hop- 
ping, while we assume a d-wave form for the BCS pairing 
hmction Ak = (Ao/2)(cosfca; — cosfcj,). The full Hamil- 
tonian for the conduction electrons reads 

= i^i3C5+7E4W<-c-W X (4) 

r 

X ($,„(r)e^K-- + <i>:„(r)e-'K-- + {x ^ y)) , 

where cr" are the Pauli spin matrices, and 7 describes the 
scattering of the Bogoliubov quasiparticles of the super- 
conductor off the SDW fluctuations. 

We computed the influence of the collective SDW 
modes on the single particle properties in the one-loop 
approximation sketched in Fig |l] with processes of sec- 
ond order in the electron-SDW coupling 7. Because of 
the breaking of translational symmetry in iSpin , the com- 
putations were numerically involved: we used lattices of 
sizes up to 32x32 with periodic boundary conditions, and 
numerically inverted matrices in real space to transform 
between the self energy and the Green's function. 

Our main results are apparent by considering the linear 
effect of the lowest order diagram which breaks trans- 
lational invariance: this is the second diagram on the 
r.h.s. of Fig. ^ and is of order Cx.yl'^- We computed 
its influence on the LDOS in a theory with JXx\i ~ 
uj'^ + (?c^ -I- A^ (J is a characteristic spin fluctuation 
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FIG. 2: The C^x,yy'^ term of as a function of the energy, 
ui, of the electron. The BCS superconductor is parameterized 
by hopping strengths t = —0.15 eV, t' = —t/4, chemical 
potential /j, — —0.135 eV which gives a doping level of about 
15%, and a size of the d-wave gap of Ao = 40 meV. Results 
are shown for different values of A, with c = 0.2 eV and ro 
on a square lattice site or bond (results for ro on a plaquette 
are similar to those on a bond). 



energy scale, and c is a spin- wave velocity) so that, after 
a Fourier transform to real space, the ^x.y propagators 
became Bessel functions of r. Following the method of 
Hoffman et al. we computed the energy (w) depen- 
dence of 



(5) 



the coefficients of the spatial Fourier transform of the 
LDOS at wavevectors K.cx,y- The Cx,yj'^ term in (H) was 
evaluated by summations over r and frequency space. Its 
values are shown in Fig ^ for a fermion band correspond- 
ing to an optimally doped cuprate. This method has the 
advantage that the values of 7, and (^x,y only modify an 
overall pre-factor of the modulation, and so do not ef- 
fect the energy dependence of the results in Fig |[ The 
change in sign of the LDOS at \uj\ < Aq arises from the 
dynamic nature of the SDW fluctuations, and the result- 
ing interference between the real and imaginary parts of 
the bare propagator and the self energy. 

We also performed computations in a different model 
in which all the diagrams in Fig |l| were summed. We 
introduced two reaZ three-component fields ipx,ya(j,T 
Re[e*^=^f ■''$x,j/a(r, t)], with the (px,ya taking momenta 
over the whole Brillouin zone of the square lattice - this 
simplifies the exact treatment of the pinning term with a 
T matrix. The bulk fluctuations were controlled by the 



susceptibilities Xi^^.„(q, w„) = 
the dispersion Wx.yi^l)'^ = 



y(q) ), where 
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FIG. 3: Grayscale plot of the quasiparticle LDOS near 
the imperfection at lattice site ro = (0,0), measured at 
uj = — O.75A0. The parameters of Hbcs are as in Fig ^ The 
SDW fluctuations are characterized by a gap of A = 5 meV, 
J = 0.15ey and a velocity of c = 0.2 eV; the SDW- 
fermion couplings is 7 = 0.1 eV and the pinning strength 
is = 0.4 eV. A similar modulation occurs for all energies 
— 2Ao < a; < — Ao/8, and significantly weaker at positive 
bias Ao/8 < < 2Ao. 




FIG. 4: As in Fig. ^ but now with a plaquette-centered defect 
at ro = (0.5,0.5). 



Ksa;,j/|) [y_7|. An important difference in this model 
is that the pinning term is now — C J (iT[(/?^^(ro, r) -f- 
(y9^^(ro, r)], with ( real; in terms of the previous formu- 
lation, this includes both |$2:,yQ,p and ^x,ya terms with 
equal modulus coefficients. This equality is a weakness 
of this second model because the physical origin of the 
_ two terms is very different, and |<I'a;,j,a(ro)|'^ is better ac- 
counted for by a shift in the local spin exciton frequency. 

Examples of the spatial dependence of the LDOS of 
this model are shown in Figs. |^ and ^. Substantial 



min(|q + Ksx,y\, |q ■ 



Fourier components at K.cx, Key were present at low en- 
ergies both for positive and negative bias. 

We also used the above results to compute (||) and the 



4 



0.01 



3 



3 

O 



-0.01 




-0.1 -0.05 0.05 0.1 

Sample bias [eV] 

FIG. 5: As in Fig ^ but for the results obtained by summing 
all the diagrams in Fig ^. The solid line is for the parameters 
in Fig. ^, while the dashed line is for 7 = 0.14 eV, 1^ = 0.2 eV. 
This result depends slightly on the system size (here 32x32 
sites) used in the numerical calculation because of the slow de- 
cay of the LDOS modulation. The peaks at ~ ±Ao — 0.04 eV 
are robust features of the theory, while the structures at other 
energies depends upon the details of the band structure and 
position of van Hove singularities. 

results are shown in Fig. |^. A robust feature of the re- 
sults in Figs ^ and ^ are the peaks in the modulus of 
the J^cx,y Fourier component near iAg. The remain- 
ing structure is influenced by both the size of the spin 
gap, A, and the coupUngs (, 7 (for a small range of ( 
values, multiple defect scattering processes can lead to 
a local low-energy bound state in the SDW propagator) , 
and depends strongly on the fermionic host band struc- 
ture, i.e., the particular form of s^. Therefore we expect 
different results for different compounds, but a detailed 
comparison with experiment is difficult, as even the bare 
LDOS seen in STM is not well represented by a two- 
dimensional tight-binding form for Sk ~ in particular the 
van-Hove singularity present in the theory is not seen 
experimentally. There are also some differences between 
the two models in Figs || and ^ at small values of the spin 
gap A: these are due to the different pinning terms near 
Tq discussed earlier, along with the modification in the 
bulk fermion spectrum associated with the first diagram 
on the r.h.s. of Fig |l| which was included only in Fig ^. 

This paper has presented a simple theory of the pin- 
ning of dynamic SDW fluctuations. Using the dynamic 
spin susceptibility (as measured in inelastic neutron scat- 
tering) and the quasiparticle dispersion as input, we pre- 
dict the energy dependence of the LDOS at the charge 
ordering wavevector, as in Figs || and ||. We hope that 
our method will ultimately allow a quantitative compar- 
ison of the spectroscopy by inelastic neutron scattering 
and by STM, and any possible disagreements should help 
expose physical effects not covered by our simple model. 

While this paper was being completed, two related 

prpnrintt; a.ppparpH 7hii pf nl (Hkp Tipf ||lfi|) His- 



reported STM measurements on Bi2Sr2CaCu208-i-a; in 
zero field, including measurements of the energy depen- 
dence of the LDOS at the wavevector {■k/2,0); the anal- 
ogous quantity for the charge order nucleated by an im- 
purity is computed here in Figs ^ and |^. 

We thank Seamus Davis for a number of discussions 
on his STM experiments and for informing us about 
his method of measuring (^. We are also grateful to 
E. Demler, A. Kapitulnik, S. Kivelson, and J. Zaanen for 
useful discussions. This research was supported by US 
NSF Grant DMR 0098226 (A.P. and S.S.) and by the 
DFG through SFB 484 (M.V.). 
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